JOURNAL OF PLANKTON RESEARCH VOLUME 33 NUMBER 1 1 PAGES 1679- 1 695 20II 



Response of Euphausia pacifica 
to small-scale shear in turbulent 
flow over a sill in a fjord 

DEBBY IANSON SUSAN E. ALLEN 2 , DAVID L. MACKAS MARK V TREVORROW 3 AND MARK C. BENFIELD 4 

'FISHERIES AND OCEANS CANADA, INSTITUTE OF OCEAN SCIENCES, PO BOX 6000, SIDNEY BC, CANADA V8L 4B2, ^DEPARTMENT OF EARTH AND OCEAN 
SCIENCES, UNIVERSITY OF BRITISH COLUMBIA, 6339 STORES RD., VANCOUVER, BC, CANADA V6T IZ4, ^DEFENCE RESEARCH AND DEVELOPMENT CANADA 
ATLANTIC, PO BOX 101 2, DARTMOUTH, NS, CANADA B2Y 3Z7 AND ^DEPARTMENT OF OCEANOGRAPHY & COASTAL SCIENCES LOUISIANA STATE UNIVERSITY, 
BATON ROUGE, LA 70803, USA 

Corresponding author: debby.ianson@dfo-mpo.gc.ca 

Received May 6, 201 1; accepted in principle July 5, 201 1; accepted for publication July 10, 201 1 
Corresponding editor: MarkJ. Gibbons 

Zooplankton in the ocean respond to visual and hydro-mechanical cues such as 
small-scale shear in turbulent flow. In addition, they form strong aggregations 
where currents intersect sloping bottoms. Strong and predictable tidal currents 
over a sill in Knight Inlet, Canada, make it an ideal location to investigate biologi- 
cal behaviour in turbulent cross-isobath flow. We examine acoustic data (38, 120 
and 200 kHz) collected there during the daylight hours, when the dominant 
zooplankters, Euphausia pacifica have descended into low light levels at ~90 m. As 
expected, these data reveal strong aggregations at the sill. However, they occur 
consistently 10— 20 m below the preferred light depth of the animals. We have 
constructed a simple model of the flow to investigate this phenomenon. Tracks of 
individual animals are traced in the flow and a variety of zooplankton behaviours 
tested. Our results indicate that the euphausiids must actively swim downward 
when they encounter the bottom boundary layer (bbl) to reproduce the observed 
downward shift in aggregation patterns. We suggest that this behaviour is cued by 
the small-scale shear in the bbl. Furthermore, this behaviour is likely to enhance 
aggregations found in strong flows at sills and on continental shelves. 

KEYWORDS: euphausiids; zooplankton behaviour; turbulence; velocity shear; 
aggregation 



INTRODUCTION 

Zooplankton in the ocean are often observed to be 
patchy spatially (Yamazaki et aL, 2003). These aggrega- 
tions occur for a variety of reasons; e.g. in response to a 
food source, as a social behaviour, to maintain position 
in convergent currents or to avoid predation (Mackas 
et aL, 1985; Genin, 2004). The presence of turbulence 



may also influence aggregation (e.g. Mackas et al, 1993; 
Franks, 2001; Tsurumi et aL, 2005). 

It has been theorized that turbulence (over spatial 
scales much larger than zooplankton) will enhance 
encounter rates of predators and prey (Rothschild and 
Osborn, 1988). Therefore, it may be expected that pre- 
dators would be attracted to regions of higher turbu- 
lence or that prey may try to avoid such regions (which 
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may also cause aggregations, e.g. Franks, 2001). It is 
also possible that animals may avoid turbulence because 
they find the velocity shear uncomfortable or because 
when the turbulence is strong enough, they may not be 
able to control their motion (e.g. Yen et al, 2008) or it 
may impact their ability to feed (e.g. Tsurumi et al., 
2005). However, it is not yet possible to make broad 
generalizations on this topic. For example, at one time 
and location, certain species of copepods have been 
observed in the turbulent layer, presumably to feed, 
while others aggregated below this layer. It was 
suggested that the animals found below may have been 
avoiding turbulence (Mackas et al, 1993). Similarly, evi- 
dence of both higher and lower feeding success in larval 
fish at elevated levels of turbulence has been found 
(Dower et al, 1988; Mackenzie et al., 1994, respectively). 

Many species of zooplankton migrate vertically to 
low light environments during the day so that predators 
cannot see them (Banse, 1964). However, zooplankton 
have been found to migrate upward to depths where 
predators are absent, in the opposite direction to what 
one would expect if they were simply avoiding light (e.g. 
Ohman, 1990). Thus, some animals alter their behav- 
iour in the presence of predators, likely due to hydro- 
mechanical (small-scale turbulent), visual or chemical 
cues (e.g. Onsrud et al, 2004). In the case of small-scale 
turbulence, many zooplankters can sense and avoid 
velocity gradients due to fluid motion caused by their 
predators (Kiorboe et al, 1999). This flight response is 
thought to be generally 'away from the turbulence', but 
not in a specific direction relative to a light gradient (i.e. 
up or down). However, Thomson and Allen (Thomson 
and Allen, 2000) observed acoustic patterns at mid- 
depths (between ~50-200m) near the continental 
shelf break, which suggested that euphausiids swam 
rapidly downward when larger predators were present 
during their daily pre-dawn descent, presumably as a 
response to one of the aforementioned cues. In this 
study, we explore the response of a specific zooplankter, 
Euphausia pacifica, to small-scale velocity shear in a 
larger-scale turbulent flow. 

Euphausia pacifica is a large zooplankter (~ 1-2 cm) 
and is generally abundant in temperate coastal Pacific 
regions (Brinton, 1962). They have large eyes that are 
extremely sensitive to light (Mauchline, 1980). They use 
their vision to maintain position within low light intensi- 
ties and undertake diel migrations on the order of 
100 m (Boden and Kampa, 1965). Euphausiids lack a 
specialized gravity-direction sensor (e.g. statocysts in 
larger decapod crustaceans) but can detect 'down' via 
light gradients, or potentially by feeling differential drag 
and rotation on their bodies if they stop swimming and 
sink in calm water (in a similar manner to the copepods 



studied by Strickler, 1982). They are sensitive to 
small-scale fluid motions (e.g. Patria and Wiese, 2004) 
and may sense small-scale shear more effectively than 
the copepods studied by Kiorboe et al. (Kiorboe et al, 

1999) because they are larger and have additional 
sensory appendages. 

Large aggregations of euphausiids have been 
observed along steep bathymetric slopes during upwel- 
ling (e.g. Simard et al, 1986; Mackas et al, 1997; Simard 
and Lavoie, 1999). Because these aggregations form and 
disperse daily, they must be caused by a rapid localized 
convergence of individuals rather than by changes in 
population size. Upwelling causes intermediate depth 
water, where these animals reside during the daylight 
hours, to be advected onshore where they encounter 
steep bathymetry and flows upward (e.g. in eastern 
boundary currents, Lentz, 1992). Simard et al. (Simard 
et al, 1986) and Simard and Mackas (Simard and 
Mackas, 1989) hypothesized that euphausiids travelling 
passively in this current will seek to maintain their pre- 
ferred light level and swim downwards when the flow 
develops a vertical component. As a result, aggregations 
form along bathymetric slopes at or slightly above the 
preferred light depth of the animals. In this research, 
we also observed aggregations of euphausiids at a 
sloping bottom, but not at their preferred light depth. 

In our study, the sloping bottom is the inner sill in 
Knight Inlet. The tidal flow over this sill is predictable 
and strong, and has been well-studied, making it an 
ideal natural laboratory. We have constructed a simple 
physical model of the upstream flow based on models 
and observation (Farmer and Armi, 1999; Cummins, 

2000) . Zooplankton were modelled as behaviour- 
capable particles embedded in this flow during the 
daytime, after the animals have migrated down. Their 
resulting paths were traced and their abundances 
sampled as a function of location. We experimented by 
implementing various animal behaviours to see which 
ones were capable of producing spatial patterns that 
reconcile with acoustic data collected at the sill. 



METHOD 

Study site 

Knight Inlet is a long (120 km), steep-sided, deep 
(>500m) fjord in mainland British Columbia, Canada 
(Fig. 1). It has strong stratification primarily due to sal- 
inity (Baker and Pond, 1995). Winds are predominately 
along the fjord but because of the stratification, their 
effect is constrained to the upper 15 m (Baker and Pond, 
1995). Our data were collected near and over the inner 
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Fig. 1 . (A) A map of Knight Inlet with the approximate location of the sill marked and (B) a longitudinal cross-section of the Knight Inlet sill, 
travelling from the seaward side, left to right at a constant latitude of 50.674°N, over a distance of 1.4 km, extracted from the echo-sounder data. 



sill (~60 m deep) roughly midway up the inlet (Fig. 1). 
The sill bottom is rocky, with rocks ranging in size from 
cobble to boulders to bedrock (Tunnicliffe and Syvistski, 
1983). Tidal amplitudes are about 3-4 m, so a large 
volume of water is exchanged diurnally (Baker and 
Pond, 1995) resulting in powerful tidal currents over the 
sill (consistendy > 60 cm s 1 near the sea-surface, 
Cummins, 2000) and often creating an internal hydraulic 
jump on the downstream side at full flow (Klymak and 
Gregg, 2001). Turbulence levels, measured as dissipation, 
6, in the strong surface flow are high, in the range of 
10" 5 - 10" 4 W kg" 1 (Klymak and Gregg, 2004). 

The sill is typical in shape relative to other sills 
formed by glaciers. It is not symmetric, rather it slopes 
more gradually on the landward (east) side of the sill 
(Fig. IB) and attains a greater depth relative to the 
seaward side. There have been times during summer 
when asymmetry in flow over the sill during the latter 
part of the tide has been observed (Klymak and Gregg, 
2004). However, the upstream flow is similar on both 
sides of the sill during the initial phase of both the ebb 
and flood tide (Klymak and Gregg, 2004). This flow is 
stratified and weakly turbulent (e = 10~ 8 -10~ 7 Wkg -1 ) 
above the bottom boundary layer (bbl). Within the bbl 
(~5-10m above the bottom), the flow is significantly 
more turbulent (e = 10 7 -10 6 Wkg consistent with 



the observations of Perlin et al. (Perlin et al, 2005) who 
show that turbulence in a shelf bbl is more than an order 
of magnitude stronger than in the flow above and bbl 
measurements made at similar current speeds in the 
nearby coastal ocean (Dewey and Crawford, 1988). This 
flow has strong horizontal velocity shear and fluid proper- 
ties are fully mixed. 

Field collection 

Intensive plankton surveys were undertaken aboard the 
CCGS Vector in November 2002 in the vicinity of the 
Knight Inlet sill (Fig. 1). The cruise was 2 weeks in 
duration and occupied several locations in the inlet 
(Trevorrow et al, 2005). Data were collected over the sill 
during 4 days (17-18 and 22-23 November). 
Two-layer estuarine circulation (Hansen and Rattray, 
1966) is low during autumn. Winds were weak to mod- 
erate during our cruises and the stratification in the 
fjord prevents their influence from penetrating below 
the upper 10-20 m (Baker and Pond, 1995). Thus, we 
were able to assume that tides dominated the physical 
circulation in our region of interest (60-100 m). 

A three-frequency (38, 120 and 200 kHz) echo- 
sounder system was used to map the spatial distributions 
of the most abundant large zooplankton, namely 
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euphausiids, amphipods and siphonophores (Trevorrow 
et al, 2005). The different frequencies were used to 
distinguish different biological target populations 
(Trevorrow et al., 2005). For example, zooplankton such 
as Euphansia pacifica exhibit an increasing acoustic target 
strength with frequency (the 38 kHz signal is much 
weaker), whereas larger animals such as fish have a 
roughly constant acoustic target strength across these 
three frequencies. Surveys were repeated across the sill 
between the 100 m isobath on the eastern side of the 
sill and ~ 1 50 m isobath on the western side over a dis- 
tance of about 2.3 km during day and night. (Figure IB 
shows the bathymetry for a typical transect.) In this 
research, the focus is on euphausiid behaviour in turbu- 
lence, so we studied the acoustic data collected during 
daylight within a few kilometres of the sill, on the 
tidal-upstream side. At night, the euphausiid scattering 
layer had migrated upward into the upper 0-20 m of 
the water column, resulting in a loss of interaction with 
the deeper flow patterns, in which turbulence is 
predictable. 

In addition to acoustic sampling, a variety of data 
were collected including a multiple net zooplankton 
sampler equipped with conductivity-temperature- 
depth and turbidity sensors (BIONESS, Bedford 
Institute of Oceanography Net Environmental 
Sampling System, Sameoto et al., 1980), and in situ 
imaging (ZOOVIS, Zooplankton Visualization System, 
Benfield et al., 2004). These data were used to 'sea- 
truth', or characterize taxonomic composition and 
target abundance in scattering layers. This sea-truthing 
was done primarily using BIONESS net tows because 
the acoustic and BIONESS sample volumes were com- 
parable (Trevorrow et al, 2005). The ZOOVIS samples 
were used only for comparison with the acoustic-derived 
abundance as its sample volume was low (about 10 4 x 
lower than either BIONESS or acoustic). Trevorrow 
et al. (Trevorrow et al., 2005) provide further detail on 



acoustic sampling methods and sea-truthing. The acous- 
tic data were then converted into individual abundances 
in spatial bins of 0.95 m (vertical) by 12 m (horizontal) 
(Appendix 1). 



Data analysis 

Every ~6— 7 h, the current over the sill changes direc- 
tion with the tide. We investigated the aggregation pat- 
terns in the acoustically derived animal abundances 
upstream of the sill on both the west and east side, 
depending on the tide. (In this study 'upstream' and 
'downstream' always refer to the tidal current.) During 
daylight hours, the euphausiid scattering layer formed a 
narrow horizontal band at 60 — 90 m depth throughout 
the Inlet (Trevorrow et al., 2005). Figure 2 shows a 
typical transect in the upstream flow, west-side in this 
case, collected at 11:19 local time at about peak flow 
(almost 3 h past slack tide). The transect provides us 
with the equivalent of a snapshot in time, showing 
animals travelling towards the sill over the slope and an 
aggregation already present at the sill. We presume that 
animals in the aggregation have been carried there 
earlier in the tide. 

In Fig. 2, the incoming scattering layer is seen at 90 m 
depth, ~ 0.5 -2 km from the sill (Fig. 2, x ~ 800- 
1200 m). The depth of the layer varied with transect, as 
expected, with isolumes changing depth depending on 
the light intensity and clarity of the water (e.g. Franks 
and Widder, 2002). We call this depth their "preferred 
light depth", 2^. Although we were unable to 
measure light levels as low as those found at ^, we 
assume that light intensities were around 
10" 4 pW cm ~ between wavelengths 475 and 480 |xm 
(Boden and Kampa, 1965; Widder and Frank, 2001). 
However, once the animals reached the sill, they consist- 
ently aggregated below ^ (Fig. 2). 
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Fig. 2. (A) A typical pattern of Euphausia pacifica abundance (individuals m ' ) near peak flow. The current is travelling from left to right 
(eastward in this case). The seafloor is dark red (corresponding to abundance greater than 300 m ), Animals aggregate spread along the seabed 
in the bottom boundary layer below their preferred light depth. (B) The same transect of euphausiid abundance with data metrics; Q, the 
incoming concentration; the preferred light depth and A^, the displacement depth of aggregation, indicated. 
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Several metrics were developed to quantitatively 
describe this aggregation and allow comparison with 
model results. The concentration factor, CF, is defined 
as the ratio of the concentration of animals in the aggre- 
gation to the average concentration of animals that are 
upstream of the sill at the preferred light depth, Q 
(from roughly 1 km upstream to the point where the sill 
depth is about 140 m, see Fig. 2B). The difference 
between the (average) preferred light depth, and the 
depth of aggregation is A£ (Fig. 2B). The fact that the 
animals occupy such a narrow vertical range (JQ) during 
the day allows the A£ to be observable. During the 
night when the animals are near the sea surface, they 
occupy a broader depth range and they may begin to 
descend slowly even when it is still dark (Thomson and 
Allen, 2000), making it much more difficult to deter- 
mine a response (such as A£) to potential turbulence. 
(There are also additional sources of backscatter near 
the sea surface, making data interpretation more 
complex, e.g. Ross et al, 2007.) 

These metrics were estimated from the abundance 
transects using a simple algorithm (see Supplementary 
Material). The vertical zone containing the majority of 
the biomass within each horizontal column of binned 
abundances was determined. The average depth, 
weighted by animal abundance, within the vertical 
band was calculated (again within each horizontal 
column of data). In the region upstream of the sill, 
these depths were averaged to yield Similarly, the 
animal weighted depths were averaged within the 
aggregation that occurred at the sill to estimate A^ (by 
difference with Supplementary Material). 

Relationships between these metrics and phase of the 
semi-diurnal tide (q>), time of day and side of sill were 
investigated to support our modelling exercise. Two 
tidal indices were used as rough indicators of the flow 
throughout the water column over the sill, the current 
index (VI) represents flow velocity and flow index (FT) 
represents the cumulative flow (from the turn of the 
current tide): 

VI = — sin <p ( 1 ) 

FT =^(l- cos <p) (2) 

where A is the amplitude of the tide (specifically the 
absolute value of the difference between tidal heights at 
the sill at slack tides straddling the sampling period, 
ranging from 2.8 to 4.5 m in our data) and <p (radians) 
is defined by the time that the aggregation was sampled 
relative to the previous slack tide. (One cycle, 2tt 
radians, is a complete tidal cycle, e.g. from one high 



tide until the next high tide.) These indices assume sinu- 
soidal variation in tidal velocity that is locally uniform 
in the vertical (barotropic) and they do not attempt to 
account for the complex baroclinic flow that may 
develop even in the upstream current at the Knight 
Inlet sill (e.g. Klymak and Gregg, 2001). All parameters 
and symbols used are listed in Appendix 2. 

Model 

We use a model to test different zooplankton behaviours 
and attempt to reproduce the patterns in the data 
described above. The model allows us to control exper- 
iments and sample zooplankton concentration at any 
time and location; something that is not possible in the 
field. It is also possible to trace the paths of individual 
animals. 

The flow in Knight Inlet is predictable and has been 
well studied. Thus, we implement a simple flow and tur- 
bulence pattern upstream of the sill (see Physical Model 
section). Individual particles (simulated zooplankters) 
are tracked in this flow using a Monte Carlo technique 
(e.g. Press et al, 1992). A large number (576 x 10 ) of 
zooplankters are followed individually in the flow field. 
The zooplankton are released steadily over 4 h (40 000 
every second). Each zooplankter enters the flow field 
over the deep plain, where the flow is not affected by 
the slope so there is no mean vertical velocity [1 km + 
25 m (random) upstream of the point where the sill 
slope begins]. They enter at a random depth between 
90 and 100 m, the observed depth of the upstream scat- 
tering layer. The motion of the zooplankters is solved 
using a fifth order Runge-Kutta scheme (Press et al., 
1992) with an adaptive time step, which in our model 
ranges from 2x10 to 10 s. The final positions are 
binned (25 m horizontal and 0.5 m vertical bins) and 
that field is then contoured. 

There are three distinct influences on a zooplankters 
path: (1) advection by the mean flow field (2) random 
motion by the turbulent flow field and (3) swimming 
motion of the zooplankter itself. Each influence is 
described in detail. 

Physical model 

The mean flow field (influence 1) is a simple, steady, 
two-dimensional representation of the modelled flow 
of Cummins (Cummins, 2000) upstream of the 
Knight Inlet sill (Fig. 3). This model was validated by 
the observations of Farmer and Armi (Farmer and 
Armi, 1999). The bbl thickness, Zb (shaded, Fig. 3) 
was specified as 10 m. The lower portion of this 
layer, the part that flows right along the seafloor, is 
termed the 'surface' layer, z s h (In this paper 'surface' 
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Fig. 3. (A) Bathymetry of model ridge on the eastward side of the sill [see equation (3)]. The mean flow field above the boundary layer is 
shown. Velocity scale vectors are shown. The vertical and horizontal scales are different. The bottom boundary layer, 1 0 m above the bottom, is 
shaded. (B) Depth profile of the mean boundary layer flow at the sill between 50 and 60 m. The bottom 1 m of the boundary layer is the 
'surface' layer. 



will always refer to the seafloor-fluid interface, rather 
than the air— sea interface.) The thickness of this 
surface layer is 10% of Zb (Stull, 1988). Details of the 
topography and the mean flow in and above the bbl 
are given in Appendix 3. All parameters and symbols 
used are listed in Appendix 2. 

Turbulence (influence 2) causes groups of animals to 
disperse. It is modelled as a diffusive process on scales 
set by the eddy viscosity, T>. In our model, each zoo- 
plankter makes a random walk in two dimensions (x and 
z) after each time step. This random walk is dependent 
on the value of T> and the time step, A^. Details of the 
values used for the eddy viscosity and the corrections to 
simple random walk required by the variation in the 
eddy viscosity (in the bbl) are given in Appendix 3. 



Zpoplankton behaviour 

We experimented with various zooplankton swimming 
behaviours (influence 3) for both east and west sill geo- 
metries (Table I). We started with the most basic; the 
null case, or no response (1) and then built in complex- 
ity until results similar to the observations were 
obtained. 

In each model scenario, the animals avoid the 
bottom, as Euphausia pacifica have been observed to in 
nature (D.L.M. and M. Tsurumi, Vancouver, using a 
remotely operated vehicle and W Hamner, Los Angeles, 
personal communication). When they come within 1 m 
of the bottom, they swim up until they are cued to 
swim down again by turbulent shear or light level 
(depending on experiment). 
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Table I: Zooplankton behaviour model 
experiments 

Model run Description 

Animals respond only to light levels, i.e. depth (exp. 1-3) 



1 . No swimming 

2. Maintain depth 

3. Maintain depth 



Passive tracer — null experiment 
At maximum swim velocity, 

zoo max = 1 0 cm s~ 1 (Genin et a/., 2005) 
Swim velocity increasing with depth away 

from z„ limited by zoo,T, ax ; 

w zoo = zootja^. - where zoOb aS ic is a 

Zfof 

comfortable euphausiid swim velocity (De 
Robertis et a/., 2003) and z xo \ is 5 m 



In addition to depth, animals respond to the bbl, i.e. small scale shear 
(exp. 4-8) 

4. Freeze and sink Euphausiids sink at 0.5 cm s~ 1 (based on 

Rubjakov, 1970) 

5. Escape response Euphausiids swim in random direction at 

zoo max 

6. Swim down Euphausiids swim downward at zoo max 



7. Swim down 



8. Swim down 



Euphausiids swim downward at zoo„ 



and 



when close to the bottom (within 5 m) they 
swim down along, and parallel to, the 
bottom 

As exp. 7 except they only are able to swim 
for 5 min at zoo max ; after which they 
become exhausted and swim at zoo basjc 



All parameters are listed in Appendix 2. 



There are not enough data available to be conclusive 
about possible patterns in A£. In general A£ is similar 
on both sides of the sill and not related to tidal indices 
[equations (1) and (2)]. As expected, CF appears to 
increase as the tide progresses and more animals are 
brought into the aggregation, i.e. with the proxy for 
integrated flow [equation (2); Table IF]. This increase is 
apparent on both sides of the sill (with the exception of 
data collected 17th November on the west side shortly 
after the turn of the tide, that has fewer animals, 
especially in the incoming scattering layer, ~5 x less 
relative to other transects). CF shows no relation to 
current index, implying that the animals are able to 
overcome the range of vertical velocities that they 
experience (Table II). 

In the absence of hydrodynamic disturbance, we 
expect the preferred depth of incoming animals to be 
related primarily to in situ light intensity. However, time 
of day is not a sufficient indicator of light intensity. 
There does not appear to be a relationship between 
<Q and time of day, or incoming velocity, in our data. 
The depth of incoming animals, ^ ranged from 60 to 
95 m, with a mean of about 85 m. In addition A£ did 
not appear to be related to time of day, or even ^j. We 
suggest that it is more likely that A£ is affected by the 
thickness of (and shear within) the bbl. 



RESULTS 

Analysis of acoustic observations 

The aggregations at the sill were consistent in their 
general pattern (Fig. 2A). On the tidal upstream side of 
the sill (whether east or west, depending on the tidal 
phase), the zooplankton scattering layer tended to dip 
10-20 m below the preferred light depth (£) as the sill 
was approached (Table II). The strongest aggregation 
was found in a dense layer oriented along and immedi- 
ately above the seabed (Fig. 2 A). Animals were concen- 
trated in these aggregations by factors ranging from 2 to 
30 (relative to C[), more often in the lower end of this 
range (Table II; mechanisms responsible for this con- 
centration are discussed with the model results). Of the 
10 transects, 3 were collected on the east side of the sill 
and 7 from the west side. The west (seaward) side of the 
sill consistently had higher concentrations of euphau- 
siids by a factor of roughly 2—4 relative to the east 
(fjord side) when comparing data collected on the same 
day. Data collected on the same tide (upstream) indicate 
that estimations of A£ were robust (varying only by 1 
and 3 m on the westward side and 6 m on the eastward 
side, Table II). The concentration factor was more vari- 
able (varying by as much as a factor of 2, Table II). 



Model experiments 

Each model experiment is discussed in the context of 
the observations and the data metrics described above. 

Passive particles (1) 

A truly passive particle, or a non-swimming zooplank- 
ter, will be carried with the flow. They rise smoothly 
from their preferred depth (90—100 m) and cross the sill 
well above (27.6-30.5 m) the bbl (Fig. 4A). A passive 
particle experiences all components of the flow. Thus, a 
group of such particles has the same divergence /conver- 
gence properties as the flow above the bbl: none. There 
is no aggregation or concentration of such particles. 

Depth maintaining behaviour (2) and (3) 
Zooplankton have been observed to swim vertically to 
maintain their preferred depth (^) = 90— 100 m in this 
case). As the flow starts over the slope and acquires a ver- 
tical velocity, the euphausiids can easily swim downward, 
against the vertical component of the flow, and maintain 
their depth (Fig. 4B). This behaviour was modelled in 
two ways, experiments (2) and (3) (Table I). The results 
of both experiments are similar. Even swimming more 
slowly (exp. 3), the animals are able to overcome the ver- 
tical component of the current and so aggregate at the 
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Table II: Table of data metrics determined from acoustic abundances (see Data Analysis) 



Date, day month 


Time, local 


Z (m) 


AZ (m) 


CF- 


if (rad) 


fa (hrs) 


VI (m) a 


FI (m) 


West side of sill 


(seaward) 
















17 November 


1 1 :24 


88 


21 


26 


0.169 


0:21 


0.34 


0.028 


18 November 


10:41 


94 


13 


30 


2.699 


5:22 


0.71 


3.1 


18 November 


10:41 


97 


16 


20 


3.041 


6:03 


0.17 


3.3 


23 November 


10:14 


94 


8.8 


2.8 


1.073 


2:03 


1.2 


0.73 


23 November 


10:57 


93 


9.2 


5.3 


1.424 


2:43 


1.4 


1.2 


23 November 


11:19 


93 


8.5 


3.5 


1.529 


2:55 


1.4 


1.3 


23 November 


11:19 


93 


8.5 


6.8 


1.623 


3:06 


1.4 


1.5 


East side of sill ( 


andward) 
















17 November 


14:52 


60 


19 


21 


1.716 


3:33 


1.9 


2.3 


1 7 November 


14:52 


58 


13 


8.6 


2.006 


3:45 


1.8 


2.8 


22 November 


14:54 


80 


20 


1.9 


0.525 


1:10 


1.1 


0.3 



The time listed is the start time of the transect, while <p and f a are calculated at the time the aggregation was sampled relative to the previous slack 
tide. (Since we only study the upstream phase of the tide ip ranges from 0 to n.) Some transects are from the same tide and at times there is more 
than one pass over the sill in the same transect (listed as separate data, but with the same start times). Two significant figures are reported for data 
metrics. 

a The velocity index is proportional to the flow velocity although it has units of length. 



sill roughly at their preferred light depth (Fig. 4B). The 
analysis yields a small A^ of -0.8 m (i.e. 0.8 m above 
their incoming depth) (Table III). The absolute error in 
our method for determining model (Supplementary 
material) is about 0.2 m. 

Aggregation mechanisms. These zooplankton are passive 
in the horizontal and thus have the same divergence/ 
convergence properties as the horizontal flow. The hori- 
zontal velocity (« a ) increases as the flow enters the 
region over the slope coming from the deep plain. 
Thus, there is a weak zooplankton divergence (they 
become more spread out in (he horizontal) over the 
slope (Fig. 4B, x ~ 2000-2500 m) followed by a strong 
convergence where the layer encounters the slope of the 
sill and zooplankton swim downwards to maintain their 
depth (Fig. 4B, x ~ 2600-2800 m). 

To quantify the aggregation processes, consider first 
an imaginary scenario with a simplified flow such that 
the mean flow is constant over the ridge (implies large 
depth) and does not slow down in the bbl but is turbu- 
lent. The zooplankton will be carried up the slope by 
the mean flow and spread through part of the bbl by 
the turbulence. The net effect is to take an initial 
biomass layer thickness, L, of say 10 m and concentrate 
it into a thinner layer, £, over the slope. This type of 
convergence, first identified over the continental slope 
(Mackas et at., 1997) and canyons (Allen et at., 2001), is 
given by L/t, where i is dependent on the properties of 
the turbulent flow in the bbl and on the behaviour of 
zooplankton when they can see the bottom. We have 
assumed that they stop trying to maintain their depth 
and swim upward once they reach the surface layer (the 
base of the bbl, just above the sea floor) so the final 



thickness is about 1.5 m. Thus, we would expect a con- 
centration factor of about 6 or 7; the simulation gives 5. 

Now consider the more realistic flow such that the 
mean flow increases over the topography and goes to 
zero against the slope. In this case, the zooplankton are 
further concentrated because they are choosing to swim 
downward across the vertical gradient of horizontal flow 
into slower moving water. The net result of the two 
effects is a concentration increase by a factor of 

iV h 

where Fj is the initial horizontal velocity of the flow and 
Vy, is the flow in the bbl (average velocity experienced 
by the zooplankton). If the incoming flow velocity is 
0.36 m s (in water of 100 m depth), then the velocity 
in the bbl at 2.0 m above bottom (at the convergence) is 
0.09 ms . Given sufficient time (>15h in our case), 
the total enhancement (or GF) in the model is 34 X . 

Over a shorter time, such as an incoming tide, the 
theoretical concentration increase for the given geome- 
try is approximately 

V,ts 

where t is the total time elapsed after the first zooplank- 
ton reach the sill and s is the slope of the ridge. Running 
the model experiment for 4 h (the time for developed 
flow of one tide in our study area) yields a CF of about 
22 on the landward side of the sill (in a 1 m thick layer) 
and 45 on the seaward side (in a 4.5 m thick layer). The 
theory above predicts a CF of 20 and 80, for landward 
and seaward geometries, respectively, if the animals have 
been flushed away from the sill during the previous tide. 
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Fig. 4. Colour contour plots of modelled concentration. The sill is shown as the white cross-hatched region. Flow goes from left to right, 
landward side of sill. (A) Passive zooplankton (exp. 1); (B) zooplankton maintain depth and swim at their maximum swim velocity (exp. 2); (C) 
zooplankton maintain depth above the boundary, in the boundary layer they freeze and passively sink (exp. 4); (D) zooplankton maintain depth 
above die boundary layer, in the boundary layer they sprint in random directions (exp. 5); (E) zooplankton maintain depth above the boundary 
layer, in the boundary layer they swim downward at maximum velocity (exp. 6); (F) zooplankton maintain depth above the boundary layer, in 
the boundary layer they swim downward and near the bottom they alter their trajectory to be downward along the topography, again at 
maximum velocity (exp. 7). 
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Table III: AZ in the modelled responses for 
various experiments 



Model run 


Depth in 


AZ 


3 depth maintain 


90.7 m 


-0.8 


4 freeze and sink 


90.6 m 


-11.8 


5 swim away in all directions 


91.0 m 


-15.3 


6 swim down 


90.9 m 


6.8 


7 swim down and along (landward) 


90.8 m 


8.5 


7 swim down and along (seaward) 


84.3 m 


13.3 



A negative AZ refers to an upward displacement. 



The assumption that the flow upslope in the bbl is negli- 
gible over the time of the aggregation is reasonable for a 
sill with a shallow slope, but is less appropriate as the 
slope becomes steeper, thus theoretical CF is overesti- 
mated on the seaward side. 

The simple depth maintenance behaviour in the 
model reproduces the observed aggregation against the 
sill and enhancements in concentration of the same 
order as in the data (see Table II, CF). However, this 
model does not result in any downward movement (A£) 
at the slope (Table III). The following experiments con- 
sider zooplankton behaviour in response to perceived 
turbulence. 



Zooplankton freeze when they encounter boundary 
layer turbulence (4) 

When model zooplankton 'freeze', or stop swimming, 
in response to turbulence, they sink at a rate of 
0.5 cm s 1 (exp. 4). This sinking rate is not enough to 
overcome the upward flow velocities on either side of 
the sill for typical tidal currents at the Knight Sill and 
so animals are swept up and over the sill (Fig. 4G). 
They would have to sink significantly faster (>2 cm s ) 
to avoid this fate. They still converge near the sill prior 
to being swept away (Fig. 4C) but there is no drop, 
rather A£ is about 12m above their preferred light 
depth (Table III). 



Zooplankton have an escape response and swim in random 
directions in response to turbulence (5) 

If zooplankton swim randomly in all directions at their 
maximum speed, they are also carried over the sill 
(exp. 5, Fig. 4D). The abundance pattern is different 
than the sinking response above. There is less aggrega- 
tion and more dispersion of the zooplankton and they 
travel faster. Similarly, they are displaced upward by 
roughly 15 m, rather than downward, at the sill 
(Table III). 



Zooplankton sprint downward in response to boundary 
layer turbulence ( 6) 

When the euphausiids escape response is to swim 
downward (as suggested by Thomson and Allen, 2000), 
the concentration patterns begin to produce a down- 
ward movement at the sill (A£ ~7 m, Table III) and 
look more similar to the observations shown in Fig. 2A 
(Fig. 4E). 

Zooplankton sprint downward and along the bottom 
in response to turbulence (7) and (8) 

When the euphausiids do not only swim vertically but 
are also allowed a horizontal component so that they 
may swim down along the bottom of the sill (Fig. 4F), 
the aggregation shifts a little further down the slope 
(A£ ~8— 13 m, Table III). We assume that since they 
are not able to sense gravity when swimming, only light 
gradients, that they choose to swim further down while 
visually avoiding the bottom (see previous section). 
They swim rapidly for a limited time from the initial 
encounter with the bbl and slow when they are close to 
the bottom (see next section). Thus, including exhaus- 
tion in the swimming behaviour (exp. 8) has only a 
small effect on the aggregation at the ridge (not shown) 
relative to exp. 7 (Fig. 4F). The centroid and deep end 
of the aggregation are nearly identical but the shallow 
end of the aggregation is slightly extended in exp. 8. Of 
all the models tested, the distribution of zooplankton in 
these experiments (on both the landward and seaward 
side of the sill; Fig. 5) are the most similar to field 
observations (Fig. 2 A, Table II). 

On the steeper seaward side, we see a larger shift 
upward because the vertical (upward) flow velocities are 
higher over the sill (Fig. 5). We used steady flow for the 
model runs described here to increase model efficiency 
and because time-varying flow did not influence the 
results. The vertical displacement obtained in time- 
varying flow depends simply on the velocity just prior to 
the time of sampling (the velocity that animals in the 
aggregation experienced as they approached the bbl). At 
peak flow, the vertical (upward) flow velocity over the sill 
is highest, causing an upward shift in the simulated 
scattering layer (as in the landward case above), so 
Aj? is slightly smaller (about 1 m relative to A£ at half 
flow). 

Path of a zooplankter 

To confirm the behaviour of individual zooplankters 
rather than the bulk properties of thousands, the paths 
of six animals are traced during the incoming tide 
for exp. 7 behaviour (Fig. 6). The zooplankton travel 
towards the sill moving up and down in the background 
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Fig. 5. Colour contour plots of modelled concentration. The sill is shown as the white cross-hatched region. Flow goes from left 
(A) On the landward side of the sill zooplankton maintain depth above the boundary layer, in the boundary layer they swim downward, 
bottom they alter their trajectory to be downward along the topography. (B) Same as (A) but for the ridge bathymetry on the western or 
side of the sill. 
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Fig. 6. Six tracks (bold traces) that illustrate depth maintenance above the boundary layer, downward swimming in the boundary layer that 
turns to parallel the slope as the slope is approached and upward swimming in the surface layer. Flow is left to right. 



turbulence swimming gently to maintain their preferred 
depth range. Where the slope begins (x = 1250 m), the 
horizontal velocity increases and they are carried 
upwards in the current, but still swim to maintain their 
depth. When they reach the bbl, they swim swiftly 
downward toward the bottom. As they approach the sill 
bottom, they slow and their tracks turn to follow the 
topography down. 



DISCUSSION 

Animal behaviour is clearly more complex than our 
simple model representations. However, our results indi- 
cate that Euphausia pacifica must actively swim down- 
ward, or to a lower light environment, when they 
encounter the bbl and small-scale velocity shear. 
Turbulence levels in the bbl in our study area are 
expected to be moderately strong, 10 7 -10 6 Wkg . 
For comparison, Euphausia pacifica themselves are able to 
generate higher turbulence when they are found in 



dense schools, 10 5 -10 4 W kg 1 (Kunze et al., 2006), 
while typical turbulence levels at depths below the wind 
mixed layer in the ocean are only about 10 9 W kg 1 . 

Because these animals are found in such a narrow ver- 
tical range during daylight hours [sensu Birch et al, 2009) 
and given the well-dehned physical flow in our study 
area, we were able to observe this behaviour in the bbl. 
They may behave similarly throughout the water 
column, for example in response to shear caused by 
wind forcing or even by a predator. The predator-escape 
response has only been observed anecdotally thus far. 
However, it seems logical that rapid responses occur at 
times because so much of a euphausiid's body is devoted 
to sprint swimming (Verity and Smetacek, 1996). On the 
other hand, it appears that Euphausia pacifica do not 
respond visually to predators in this manner, at least 
when the fish are more than 20 cm away (De Robertis 
et al, 2003). A downward response to turbulence caused 
by wind forcing in the upper mixed layer would be diffi- 
cult to observe, but may be of consequence to ecological 
models of grazing by zooplankton. 
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Strong aggregations of zooplankton have been found 
at other sea slopes during the day where currents flow 
upwards and the animals are able to maintain their light 
depth (Genin, 2004). In eastern boundary currents, the 
cross-isobath, or upwelling, flow is relatively weak. The 
strong flow is in the along-isobath (along-shore) direction 
(Lentz, 1992). Given the weak flow, the documented 
aggregations in these regions (e.g. Simard and Mackas, 
1989) are more dense than might be expected for diel 
migrators. However, the currents in the along-shore 
direction are strong and consistent enough to produce 
boundary layer shear even if the on-shore current is not. 
Our results suggest that the aggregations occurring at 
these locations may be further concentrated, in-line with 
observations, by the downward swimming behaviour 
that we describe, because travelling into slower moving 
water causes an additional convergence. 

The euphausiids do not appear to be moving down- 
ward in response to predators. At times larger animals, 
presumably fish, are found in the bbl at the sill in dense 
aggregations (seen in the lower frequency, 38 kHz, 
data). However, these large scatterers were usually well 
below the euphausiids' preferred light depth (that 
ranged from 60 to 100 m) at around 140 m. 

The motivation for the model study relies on the 
acoustic data. The conversion from acoustic backscatter 
to animal abundance, and especially to a specific zoo- 
plankter, is complex and requires many assumptions (e.g. 
Benfield et al, 1998). We were fortunate to have high 
abundances in 2002 and in situ net samples to satisfy 
some of the important assumptions. For example, the net 
samples showed that a single species dominated the 
intense scattering layer (Fig. 2) and that its size distri- 
bution was Gaussian (Fig. 3, Trevorrow et al., 2005). The 
38 kHz data confirmed that the contribution of larger 
animals such as fish was minimal within the scattering 
layer and the abundance algorithm was designed to filter 
out that portion of the signal (Trevorrow et al., 2005). In 
addition, the ZOOVIS data suggest that our acoustic 
abundances are underestimated by roughly a factor of 2 
which could be a result of changes in swimming orien- 
tation (Trevorrow et al, 2005). Thus, we caution the 
reader that our quantitative estimates from data should 
be allowed a margin of error. The primary function of 
this analysis (Table II) is to allow comparison between 
the data and model results. 

The observations show the higher concentrations of 
animals at the sill that both physical theory and model 
predict. On the landward side, all three measures of 
this enhancement (CF) agree, at the high end of the 
range in the data (see Table II, and Model results). On 
the seaward side, model and theory over-predict 
CF (see Model results) relative to the data. The model is 



a factor of 1.5 x higher than the high end of the data 
range (30, Table II) while the theoretical CF is even 
higher, almost double that of the model. We never 
observed CF higher than 30 x in the data and consider 
that model and theory may allow overcrowding that 
does not occur in nature. The model also has the 
advantage of injecting many animals at consistent inter- 
vals so that clear and repeatable responses are possible. 
In the observations, we do not know how many animals 
have been swept toward the sill prior to the snap-shot 
that the transect provides. 

Exploration of relationships between these metrics 
(CF and A£) with current speed and integrated flow [or 
time in tidal cycle, see equations (1) and (2)] are incon- 
clusive because there are not enough data. However, 
they do suggest that CF increases as the tidal cycle pro- 
gresses and more animals are carried toward the sill as 
one would expect. The bbl may become thicker and 
more turbulent as the mean flow increases (lagging the 
mean flow) potentially influencing both A£ and CF. 
Model runs show that decreasing bbl thickness from 
10 m, used in our standard model runs, to 5 m (after 
Cummins, 2000) led to a similar (approximately 50%) 
decrease in the downward shift of the aggregation, 
while CF was similar. Thus we might expect to 
increase as the tidal cycle progresses. 

We assume that the depth of isolumes remains con- 
stant in the region upstream of the sill. The fresher 
surface layer, which may absorb more light, becomes 
significantly thicker over the top of the sill. However, it 
is constant upstream of the sill, deepening slightly by 
~10% in the region of aggregation (Klymak and 
Gregg, 2004). Regardless, if there were more light 
absorption above the aggregation, relative to the region 
upstream, then we would expect the euphausiids to be 
found higher in the water column to maintain their 
light depth, rather than lower, as the sill is approached. 

The data are consistent on either side of the sill 
despite differences in slope and possibly in flow 
characteristics during the later part of the tide. On the 
west side, water around the preferred light depth ~ 
90 m) may slow as the tide progresses due to a dense 
(saltier than on the landward side) body of water often 
present there (observed during summer by Klymak 
and Gregg, 2003). More data were collected on the 
west side of the sill, and two of those transects (col- 
lected on 11/18, Table II) occur when the water at 
^ at the sill may be moving slowly. In both cases, 
aggregations are seen below If the water is no 
longer turbulent, we do not know why Euphausia pacifica 
remain below their preferred light depth, unless there 
is a lag time before they venture back up after experi- 
encing hydrodynamical disturbances. 
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CONCLUSIONS 

We propose a new mechanism to enhance zooplankton 
aggregation and a downward swimming response to 
small-scale shear by zooplankton. Bio-physical aggrega- 
tion in the marine environment in response to upwel- 
ling flow over topographical features has been observed 
and documented (Genin, 2004). Previous studies indi- 
cate that during the day, light-sensitive animals choose 
to maintain their light depth (Genin, 2004). However, 
in the strong tidal flows at the Knight Inlet sill Euphausia 
pacifica aggregate below this depth. We suggest that the 
downward shift in aggregation is caused by a behaviour- 
al response to small-scale velocity shear in the turbulent 
bbl. This response causes the concentration of the 
animals to be further enhanced as they move into even 
slower moving water. In addition, our model shows that 
animals must actively swim down to a lower light 
environment to produce a downward shift in concen- 
tration patterns in such a strong current. It is not poss- 
ible to produce such a pattern by freezing and sinking 
or by swimming fast in random directions. We suggest 
that this response by euphausiids is not confined to our 
study area and that it is likely an evolved response to 
small-scale velocity shear (turbulence), caused by wind, 
tides or predators, throughout the water column. 

SUPPLEMENTARY DATA 

Supplementary data can be found online at http:// 
plankt. oxfordj ournals .org. 
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APPENDIX: 1: ACOUSTIC DATA 
PROCESSING 

The raw echo-sounder data were averaged vertically 
into 0.95 m depth bins, and horizontally by 4 pings 
(4 s or about 12 m in horizontal distance at a speed 
of 6 knots). These data were first converted to volu- 
metric backscatter strength and the standard correc- 
tions made (Trevorrow et al., 2005). The calibration 
data required for the corrections were collected while 
the ship was anchored in Hoeya Sound (beside the 
sill, 50°42.00'N, 125°58.00'\y Fig. 1A) and are esti- 
mated to be accurate to + 0.8 dB, or a relative error 
of 20%. In our data, the observed backscatter within 
the scattering layer varied by a factor of roughly 

2- 10, while outside of this layer the backscatter was 

3— 4 orders of magnitude less. 

The BIONESS and 38 kHz acoustic data indicated 
that the mid-depth (90 m) sound scattering layer near 
the sill was strongly dominated by Euphausia pacifica 
(Trevorrow et al., 2005). The size and age distribution 
of these E. pacifica was unimodal, dominated by young 
adults born the previous spring. Elsewhere in the inlet 
there were times and places where most of the 
biomass in the deep scattering layer was made up of 
other zooplankters, such as siphonophores or large 
copepods, and fish (Trevorrow et al., 2005). Scattering 
in the upper layer near the sill is caused primarily by 
physical microstructure (Ross and Lueck, 2003). Thus, 
in situ abundances were estimated by dividing the 
volumetric backscatter cross-section estimates by the 
population-averaged backscatter cross-section for an 
individual euphausiid at each of the echo-sounder 
frequencies. 

The population averaged acoustic scattering model 
from Stanton and Chu (Stanton and Chu, 2000) was 
the most appropriate to our study (Trevorrow et al., 
2005). Euphausiids are represented by bent fluid cylin- 
ders. Distributions of their length and orientation are 
explicitly averaged within the model. We used the 
BIONESS data to determine the distribution of cylin- 
der dimensions; length = (15.9 + 1.3) mm (mean, 
standard deviation) with a length to radius ratio of 15 
(M. Galbraith, unpublished data). We assume that the 
radius of curvature of the cylinders is 3 x the body 



length and that their orientation is horizontal (0°) + 
30° (mean, standard deviation). The sound speed and 
density contrast ratios (zooplankter vs. seawater) were 
assumed to be 1.025 and 1.045 (following Foote et al., 
1990). Thus, the model predicts average backscatter 
target strengths per animal of —97.8, —83.1 and 
-78.9 dB (re 1 m 2 ) for 38, 120 and 200 kHz, 
respectively. Abundances predicted by the model were 
quantitatively verified with simultaneous BIONESS 
trawls (Trevorrow et al., 2005). 



APPENDIX 2: MODEL 
PARAMETERS AND DATA 
METRICS 

All parameters and symbols used are listed in Table IV 

Table IV: Model parameters and data metrics 
defined with units 

Physical model parameters 
V - Eddy viscosity 

h 0 m Average water column depth upstream of sill, 

constant 

h b m Water column depth over sill, variable 



u 
\t 

u 
u 

u B 

w m 

X 

z 
A> 

Zooplankton behaviour parameters 



W zoo 


m s _1 


Vertical component of zooplankter swim velocity 


Zj 


m 


Preferred light depth of zooplankter 




m 


Zooplankter's "acceptable" distance away from z. 


z00 basic 


m s _1 


Comfortable zooplankter swimming speed 


zoo max 


m s _1 


Maximum zooplankter swimming speed 


Data metrics 




<P 


rad 


Phase of tide 


CF 




Concentration factor 


Q 


# rrT 3 


Concentration of animals approaching the sill 


11 


m 


Flow index 


t a 


h 


Time aggregation was sampled relative to previous 






slack tide 


VI 


m 


Velocity index 


AZ 


m 


Displacement depth of aggregation below Z, 


Z| 


m 


Preferred light depth of animals approaching 



the sill 



Eddy viscosity 
m Average water column depth upstream of sill, 

constant 

m Water column depth over sill, variable 

m Water column depth at top of sill, constant 

m 2 s~ 1 Eddy viscosity in interior flow 

Von Karman constant 
m Horizontal location of beginning of sill slope 

m Horizontal location of top of sill 

s Time step in model 

m s _1 Mean flow velocity upstream of sill, constant 
m s~ 1 General mean horizontal velocity 
m s _1 Horizontal velocity in the mean flow over sill, 

variable 
m s _1 Friction velocity 
m s~ 1 Mean vertical velocity 
m General horizontal coordinate 

m General vertical coordinate 

m Thickness of bottom boundary layer (bbl) 

m Thickness of 'surface' layer, bottom of bbl 
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APPENDIX 3: PHYSICAL MODEL 
Topography 

The model replicates the main features of the topogra- 
phy on each side of the sill. In general, there is a deep 
plain, followed by a linear slope and a shallow plain 
(Fig. 3A). On the east (landward) side, the bottom 
depth, hb, is given by 



ho 



hbE(x) = < K - (k) - /'min) 



fin 



x — L\ 
L> - L 



x < L\ 
L x <x< L 2 
x > Li 



(3) 



where /z mm = 60 m, h Q = 200 m, L x = 1000 m, L 2 = 
3250 m and x is the horizontal location increasing as 
the sill is approached (Fig. 3A). On the west (seaward) 
side, the slope has two sections and is steeper, especially 
near the top of the sill (Fig. IB). The basin is not quite 
as deep. The bottom depth is given by 



Abw(*) 



/'ii 

ha - (h 0 - Ih) 



hi ~ (hi - hmin) 



x ■ 



Li 



x < Li 

Li < X < Ly, 
Lb < X < Li 

x > Li 



(4) 



where again /z min = 60 m, while h Q = 150 m, hi - 
130 m, Li = 2000 m, L b = 2500 m and L 2 = 2880 m. 



Mean flow field 

The mean flow is based on the Cummins (Cummins, 
2000) model. It is steady and conserves volume as it 
passes over the slope. The horizontal velocity w a in the 
model above the bottom layer (Fig. 3A) is: 



u&(x) = U 



hp 
hb(x) 



(3) 



where t/=0.18ms 1 is flow speed approaching the 
sill on the east side, and is about 10% higher on the 
west side (because the water column on the west side is 
shallower, so the same volume of water flows more 
quickly). 

Friction velocity, u r , is a measure of the surface stress. 
It is based on the velocity and eddy viscosity (k m = 
0.01 m 2 s ) above the bbl (bottom 10 m, shaded 
Fig. 3A) in the mean flow (Cummins, 2000) so that u* = 
(kmU^x) I ' Zb) The flow velocity, u(x, z), in the surface 



layer (bottom 1 0% of bbl) was assumed to vary linearly 
with depth (Fig. 3B) for simplicity (using a conventional 
log layer makes no difference to the zooplankton advec- 
tion/ aggregation patterns); 

u(x, z) = u*(x) Z+ b ^ , z < -h h (x) + Zd (6) 

where z is the vertical coordinate measured upward 
from the air— sea interface (so that z is always negative 
in the ocean). Above the surface layer in the bbl, the 
flow is quadratic and the velocities at the boundaries of 
adjacent layers are continuous. 



h h (x) 



l(x,z) = K a (*) - («a(*) - "*(#)) 

£b - Z s l 

- h h (x) + Zs i<Z< -K(x) + Zb (7) 
and above the bbl 

u(x, z) = -K(x) + Zb < Z (8) 

The mean vertical velocity w m , set by conservation of 
mass, is: 



0 



x < L\ , x > Li 



^ *) = < -„(,, Z ) h -f^^- Li<x<L 2 
L-2 - L\ h b (x) 



Turbulent flow 

In the free flow, the eddy viscosity is taken as a 
constant. 



V = k mi —hb + Zb < Z 



(10) 



The eddy viscosity in the bbl was modelled using the 
O'Brien (O'Brien, 1970) formulation: 



-hb+Zb-Z 
Zb — Zd 



- hb + zj < z< ~H h + Zb 



KU*Z s \-kn 

Zb ^ Zsl 



(11) 



where k is the dimensionless von Karman constant 
(0.4). In the surface layer, the eddy viscosity is given as: 

V= Kiw( Z + h h ),z < -h + ZA (12) 

The maximum turbulent velocity shear occurs at about 
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40% of the depth of the bbl. The random walk is 
denned as a vertical step and a horizontal step of a 
random number (R, mean 0 and variance 1) multiplied 
by (2VAt) 1/2 . In the bbl, the random walk must be cor- 
rected because of the non-constant diffusivity (Hunter 
et al, 1993; Visser, 1997). Thus, the turbulent velocity 
was taken as 











f R 






M. 



where 2?' is the derivative of the diffusivity in the ver- 
tical or horizontal direction for the vertical or hori- 
zontal step, respectively. Using an adaptive step-size 
can be an issue for non-uniform diffusivity (Brickman 
and Smith, 2002) but was not found to be here. This 
formulation for the turbulent velocities for the 
Runge— Kutta scheme maintained an initially well- 
mixed particle field (as required — see Thomson, 
1984). The coefficient on the derivative of diffusivity 
is, on average, half that needed for a simple Euler 
scheme (Visser, 1997). 
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